3.4.76 \(\int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx\) [376]

3.4.76.1 Optimal result
3.4.76.2 Mathematica [C] (warning: unable to verify)
3.4.76.3 Rubi [A] (verified)
3.4.76.4 Maple [C] (verified)
3.4.76.5 Fricas [B] (verification not implemented)
3.4.76.6 Sympy [F]
3.4.76.7 Maxima [F]
3.4.76.8 Giac [B] (verification not implemented)
3.4.76.9 Mupad [F(-1)]

3.4.76.1 Optimal result

Integrand size = 25, antiderivative size = 97 \[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\frac {(a-2 b) \arctan \left (\frac {\sqrt {a-b} \sinh (e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}}\right )}{2 (a-b)^{3/2} f}+\frac {\text {sech}(e+f x) \sqrt {a+b \sinh ^2(e+f x)} \tanh (e+f x)}{2 (a-b) f} \]

output
1/2*(a-2*b)*arctan(sinh(f*x+e)*(a-b)^(1/2)/(a+b*sinh(f*x+e)^2)^(1/2))/(a-b 
)^(3/2)/f+1/2*sech(f*x+e)*(a+b*sinh(f*x+e)^2)^(1/2)*tanh(f*x+e)/(a-b)/f
 
3.4.76.2 Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 5 vs. order 3 in optimal.

Time = 6.25 (sec) , antiderivative size = 417, normalized size of antiderivative = 4.30 \[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\frac {\text {sech}^3(e+f x) \left (1+\frac {b \sinh ^2(e+f x)}{a}\right ) \tanh (e+f x) \left (45 a \arcsin \left (\sqrt {\frac {(a-b) \tanh ^2(e+f x)}{a}}\right )+30 b \arcsin \left (\sqrt {\frac {(a-b) \tanh ^2(e+f x)}{a}}\right ) \sinh ^2(e+f x)+16 a \operatorname {Hypergeometric2F1}\left (2,3,\frac {7}{2},\frac {(a-b) \tanh ^2(e+f x)}{a}\right ) \sqrt {\frac {\text {sech}^2(e+f x) \left (a+b \sinh ^2(e+f x)\right )}{a}} \left (\frac {(a-b) \tanh ^2(e+f x)}{a}\right )^{5/2}+16 b \operatorname {Hypergeometric2F1}\left (2,3,\frac {7}{2},\frac {(a-b) \tanh ^2(e+f x)}{a}\right ) \sinh ^2(e+f x) \sqrt {\frac {\text {sech}^2(e+f x) \left (a+b \sinh ^2(e+f x)\right )}{a}} \left (\frac {(a-b) \tanh ^2(e+f x)}{a}\right )^{5/2}-45 a \sqrt {\frac {(a-b) \text {sech}^2(e+f x) \left (a+b \sinh ^2(e+f x)\right ) \tanh ^2(e+f x)}{a^2}}-30 b \sinh ^2(e+f x) \sqrt {\frac {(a-b) \text {sech}^2(e+f x) \left (a+b \sinh ^2(e+f x)\right ) \tanh ^2(e+f x)}{a^2}}\right )}{30 a f \sqrt {a+b \sinh ^2(e+f x)} \sqrt {\frac {\text {sech}^2(e+f x) \left (a+b \sinh ^2(e+f x)\right )}{a}} \left (\frac {(a-b) \tanh ^2(e+f x)}{a}\right )^{3/2}} \]

input
Integrate[Sech[e + f*x]^3/Sqrt[a + b*Sinh[e + f*x]^2],x]
 
output
(Sech[e + f*x]^3*(1 + (b*Sinh[e + f*x]^2)/a)*Tanh[e + f*x]*(45*a*ArcSin[Sq 
rt[((a - b)*Tanh[e + f*x]^2)/a]] + 30*b*ArcSin[Sqrt[((a - b)*Tanh[e + f*x] 
^2)/a]]*Sinh[e + f*x]^2 + 16*a*Hypergeometric2F1[2, 3, 7/2, ((a - b)*Tanh[ 
e + f*x]^2)/a]*Sqrt[(Sech[e + f*x]^2*(a + b*Sinh[e + f*x]^2))/a]*(((a - b) 
*Tanh[e + f*x]^2)/a)^(5/2) + 16*b*Hypergeometric2F1[2, 3, 7/2, ((a - b)*Ta 
nh[e + f*x]^2)/a]*Sinh[e + f*x]^2*Sqrt[(Sech[e + f*x]^2*(a + b*Sinh[e + f* 
x]^2))/a]*(((a - b)*Tanh[e + f*x]^2)/a)^(5/2) - 45*a*Sqrt[((a - b)*Sech[e 
+ f*x]^2*(a + b*Sinh[e + f*x]^2)*Tanh[e + f*x]^2)/a^2] - 30*b*Sinh[e + f*x 
]^2*Sqrt[((a - b)*Sech[e + f*x]^2*(a + b*Sinh[e + f*x]^2)*Tanh[e + f*x]^2) 
/a^2]))/(30*a*f*Sqrt[a + b*Sinh[e + f*x]^2]*Sqrt[(Sech[e + f*x]^2*(a + b*S 
inh[e + f*x]^2))/a]*(((a - b)*Tanh[e + f*x]^2)/a)^(3/2))
 
3.4.76.3 Rubi [A] (verified)

Time = 0.30 (sec) , antiderivative size = 101, normalized size of antiderivative = 1.04, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.200, Rules used = {3042, 3669, 296, 291, 216}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {1}{\cos (i e+i f x)^3 \sqrt {a-b \sin (i e+i f x)^2}}dx\)

\(\Big \downarrow \) 3669

\(\displaystyle \frac {\int \frac {1}{\left (\sinh ^2(e+f x)+1\right )^2 \sqrt {b \sinh ^2(e+f x)+a}}d\sinh (e+f x)}{f}\)

\(\Big \downarrow \) 296

\(\displaystyle \frac {\frac {(a-2 b) \int \frac {1}{\left (\sinh ^2(e+f x)+1\right ) \sqrt {b \sinh ^2(e+f x)+a}}d\sinh (e+f x)}{2 (a-b)}+\frac {\sinh (e+f x) \sqrt {a+b \sinh ^2(e+f x)}}{2 (a-b) \left (\sinh ^2(e+f x)+1\right )}}{f}\)

\(\Big \downarrow \) 291

\(\displaystyle \frac {\frac {(a-2 b) \int \frac {1}{1-\frac {(b-a) \sinh ^2(e+f x)}{b \sinh ^2(e+f x)+a}}d\frac {\sinh (e+f x)}{\sqrt {b \sinh ^2(e+f x)+a}}}{2 (a-b)}+\frac {\sinh (e+f x) \sqrt {a+b \sinh ^2(e+f x)}}{2 (a-b) \left (\sinh ^2(e+f x)+1\right )}}{f}\)

\(\Big \downarrow \) 216

\(\displaystyle \frac {\frac {(a-2 b) \arctan \left (\frac {\sqrt {a-b} \sinh (e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}}\right )}{2 (a-b)^{3/2}}+\frac {\sinh (e+f x) \sqrt {a+b \sinh ^2(e+f x)}}{2 (a-b) \left (\sinh ^2(e+f x)+1\right )}}{f}\)

input
Int[Sech[e + f*x]^3/Sqrt[a + b*Sinh[e + f*x]^2],x]
 
output
(((a - 2*b)*ArcTan[(Sqrt[a - b]*Sinh[e + f*x])/Sqrt[a + b*Sinh[e + f*x]^2] 
])/(2*(a - b)^(3/2)) + (Sinh[e + f*x]*Sqrt[a + b*Sinh[e + f*x]^2])/(2*(a - 
 b)*(1 + Sinh[e + f*x]^2)))/f
 

3.4.76.3.1 Defintions of rubi rules used

rule 216
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*A 
rcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a 
, 0] || GtQ[b, 0])
 

rule 291
Int[1/(Sqrt[(a_) + (b_.)*(x_)^2]*((c_) + (d_.)*(x_)^2)), x_Symbol] :> Subst 
[Int[1/(c - (b*c - a*d)*x^2), x], x, x/Sqrt[a + b*x^2]] /; FreeQ[{a, b, c, 
d}, x] && NeQ[b*c - a*d, 0]
 

rule 296
Int[((a_) + (b_.)*(x_)^2)^(p_)*((c_) + (d_.)*(x_)^2)^(q_), x_Symbol] :> Sim 
p[(-b)*x*(a + b*x^2)^(p + 1)*((c + d*x^2)^(q + 1)/(2*a*(p + 1)*(b*c - a*d)) 
), x] + Simp[(b*c + 2*(p + 1)*(b*c - a*d))/(2*a*(p + 1)*(b*c - a*d))   Int[ 
(a + b*x^2)^(p + 1)*(c + d*x^2)^q, x], x] /; FreeQ[{a, b, c, d, q}, x] && N 
eQ[b*c - a*d, 0] && EqQ[2*(p + q + 2) + 1, 0] && (LtQ[p, -1] ||  !LtQ[q, -1 
]) && NeQ[p, -1]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3669
Int[cos[(e_.) + (f_.)*(x_)]^(m_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^ 
(p_.), x_Symbol] :> With[{ff = FreeFactors[Sin[e + f*x], x]}, Simp[ff/f   S 
ubst[Int[(1 - ff^2*x^2)^((m - 1)/2)*(a + b*ff^2*x^2)^p, x], x, Sin[e + f*x] 
/ff], x]] /; FreeQ[{a, b, e, f, p}, x] && IntegerQ[(m - 1)/2]
 
3.4.76.4 Maple [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 1.52 (sec) , antiderivative size = 35, normalized size of antiderivative = 0.36

method result size
default \(\frac {\operatorname {`\,int/indef0`\,}\left (\frac {1}{\cosh \left (f x +e \right )^{4} \sqrt {a +b \sinh \left (f x +e \right )^{2}}}, \sinh \left (f x +e \right )\right )}{f}\) \(35\)
risch \(\text {Expression too large to display}\) \(108002\)

input
int(sech(f*x+e)^3/(a+b*sinh(f*x+e)^2)^(1/2),x,method=_RETURNVERBOSE)
 
output
`int/indef0`(1/cosh(f*x+e)^4/(a+b*sinh(f*x+e)^2)^(1/2),sinh(f*x+e))/f
 
3.4.76.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 693 vs. \(2 (85) = 170\).

Time = 0.32 (sec) , antiderivative size = 1503, normalized size of antiderivative = 15.49 \[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\text {Too large to display} \]

input
integrate(sech(f*x+e)^3/(a+b*sinh(f*x+e)^2)^(1/2),x, algorithm="fricas")
 
output
[-1/4*(((a - 2*b)*cosh(f*x + e)^4 + 4*(a - 2*b)*cosh(f*x + e)*sinh(f*x + e 
)^3 + (a - 2*b)*sinh(f*x + e)^4 + 2*(a - 2*b)*cosh(f*x + e)^2 + 2*(3*(a - 
2*b)*cosh(f*x + e)^2 + a - 2*b)*sinh(f*x + e)^2 + 4*((a - 2*b)*cosh(f*x + 
e)^3 + (a - 2*b)*cosh(f*x + e))*sinh(f*x + e) + a - 2*b)*sqrt(-a + b)*log( 
((a - 2*b)*cosh(f*x + e)^4 + 4*(a - 2*b)*cosh(f*x + e)*sinh(f*x + e)^3 + ( 
a - 2*b)*sinh(f*x + e)^4 - 2*(3*a - 2*b)*cosh(f*x + e)^2 + 2*(3*(a - 2*b)* 
cosh(f*x + e)^2 - 3*a + 2*b)*sinh(f*x + e)^2 - 2*sqrt(2)*(cosh(f*x + e)^2 
+ 2*cosh(f*x + e)*sinh(f*x + e) + sinh(f*x + e)^2 - 1)*sqrt(-a + b)*sqrt(( 
b*cosh(f*x + e)^2 + b*sinh(f*x + e)^2 + 2*a - b)/(cosh(f*x + e)^2 - 2*cosh 
(f*x + e)*sinh(f*x + e) + sinh(f*x + e)^2)) + 4*((a - 2*b)*cosh(f*x + e)^3 
 - (3*a - 2*b)*cosh(f*x + e))*sinh(f*x + e) + a - 2*b)/(cosh(f*x + e)^4 + 
4*cosh(f*x + e)*sinh(f*x + e)^3 + sinh(f*x + e)^4 + 2*(3*cosh(f*x + e)^2 + 
 1)*sinh(f*x + e)^2 + 2*cosh(f*x + e)^2 + 4*(cosh(f*x + e)^3 + cosh(f*x + 
e))*sinh(f*x + e) + 1)) - 2*sqrt(2)*((a - b)*cosh(f*x + e)^2 + 2*(a - b)*c 
osh(f*x + e)*sinh(f*x + e) + (a - b)*sinh(f*x + e)^2 - a + b)*sqrt((b*cosh 
(f*x + e)^2 + b*sinh(f*x + e)^2 + 2*a - b)/(cosh(f*x + e)^2 - 2*cosh(f*x + 
 e)*sinh(f*x + e) + sinh(f*x + e)^2)))/((a^2 - 2*a*b + b^2)*f*cosh(f*x + e 
)^4 + 4*(a^2 - 2*a*b + b^2)*f*cosh(f*x + e)*sinh(f*x + e)^3 + (a^2 - 2*a*b 
 + b^2)*f*sinh(f*x + e)^4 + 2*(a^2 - 2*a*b + b^2)*f*cosh(f*x + e)^2 + 2*(3 
*(a^2 - 2*a*b + b^2)*f*cosh(f*x + e)^2 + (a^2 - 2*a*b + b^2)*f)*sinh(f*...
 
3.4.76.6 Sympy [F]

\[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\int \frac {\operatorname {sech}^{3}{\left (e + f x \right )}}{\sqrt {a + b \sinh ^{2}{\left (e + f x \right )}}}\, dx \]

input
integrate(sech(f*x+e)**3/(a+b*sinh(f*x+e)**2)**(1/2),x)
 
output
Integral(sech(e + f*x)**3/sqrt(a + b*sinh(e + f*x)**2), x)
 
3.4.76.7 Maxima [F]

\[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\int { \frac {\operatorname {sech}\left (f x + e\right )^{3}}{\sqrt {b \sinh \left (f x + e\right )^{2} + a}} \,d x } \]

input
integrate(sech(f*x+e)^3/(a+b*sinh(f*x+e)^2)^(1/2),x, algorithm="maxima")
 
output
integrate(sech(f*x + e)^3/sqrt(b*sinh(f*x + e)^2 + a), x)
 
3.4.76.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 693 vs. \(2 (85) = 170\).

Time = 0.43 (sec) , antiderivative size = 693, normalized size of antiderivative = 7.14 \[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\frac {{\left (\frac {{\left (a - 2 \, b\right )} \arctan \left (-\frac {\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b} + \sqrt {b}}{2 \, \sqrt {a - b}}\right )}{{\left (a e^{\left (4 \, e\right )} - b e^{\left (4 \, e\right )}\right )} \sqrt {a - b}} - \frac {2 \, {\left ({\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )}^{3} a - 2 \, {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )}^{3} b - 5 \, {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )}^{2} a \sqrt {b} + 2 \, {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )}^{2} b^{\frac {3}{2}} - 4 \, {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )} a^{2} - {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )} a b + 2 \, {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )} b^{2} - 4 \, a^{2} \sqrt {b} + 5 \, a b^{\frac {3}{2}} - 2 \, b^{\frac {5}{2}}\right )}}{{\left ({\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )}^{2} + 2 \, {\left (\sqrt {b} e^{\left (2 \, f x + 2 \, e\right )} - \sqrt {b e^{\left (4 \, f x + 4 \, e\right )} + 4 \, a e^{\left (2 \, f x + 2 \, e\right )} - 2 \, b e^{\left (2 \, f x + 2 \, e\right )} + b}\right )} \sqrt {b} + 4 \, a - 3 \, b\right )}^{2} {\left (a e^{\left (4 \, e\right )} - b e^{\left (4 \, e\right )}\right )}}\right )} e^{\left (4 \, e\right )}}{f} \]

input
integrate(sech(f*x+e)^3/(a+b*sinh(f*x+e)^2)^(1/2),x, algorithm="giac")
 
output
((a - 2*b)*arctan(-1/2*(sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4*f*x + 4*e) + 
 4*a*e^(2*f*x + 2*e) - 2*b*e^(2*f*x + 2*e) + b) + sqrt(b))/sqrt(a - b))/(( 
a*e^(4*e) - b*e^(4*e))*sqrt(a - b)) - 2*((sqrt(b)*e^(2*f*x + 2*e) - sqrt(b 
*e^(4*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 2*b*e^(2*f*x + 2*e) + b))^3*a - 2 
*(sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 
 2*b*e^(2*f*x + 2*e) + b))^3*b - 5*(sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4* 
f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 2*b*e^(2*f*x + 2*e) + b))^2*a*sqrt(b) + 
 2*(sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) 
 - 2*b*e^(2*f*x + 2*e) + b))^2*b^(3/2) - 4*(sqrt(b)*e^(2*f*x + 2*e) - sqrt 
(b*e^(4*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 2*b*e^(2*f*x + 2*e) + b))*a^2 - 
 (sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 
 2*b*e^(2*f*x + 2*e) + b))*a*b + 2*(sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4* 
f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 2*b*e^(2*f*x + 2*e) + b))*b^2 - 4*a^2*s 
qrt(b) + 5*a*b^(3/2) - 2*b^(5/2))/(((sqrt(b)*e^(2*f*x + 2*e) - sqrt(b*e^(4 
*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 2*b*e^(2*f*x + 2*e) + b))^2 + 2*(sqrt( 
b)*e^(2*f*x + 2*e) - sqrt(b*e^(4*f*x + 4*e) + 4*a*e^(2*f*x + 2*e) - 2*b*e^ 
(2*f*x + 2*e) + b))*sqrt(b) + 4*a - 3*b)^2*(a*e^(4*e) - b*e^(4*e))))*e^(4* 
e)/f
 
3.4.76.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\text {sech}^3(e+f x)}{\sqrt {a+b \sinh ^2(e+f x)}} \, dx=\int \frac {1}{{\mathrm {cosh}\left (e+f\,x\right )}^3\,\sqrt {b\,{\mathrm {sinh}\left (e+f\,x\right )}^2+a}} \,d x \]

input
int(1/(cosh(e + f*x)^3*(a + b*sinh(e + f*x)^2)^(1/2)),x)
 
output
int(1/(cosh(e + f*x)^3*(a + b*sinh(e + f*x)^2)^(1/2)), x)